Double dark state cooling in a three level system 
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A detailed study of a robust and fast laser cooling scheme on a three level system is presented. 
A special laser configuration, applicable to trapped ions, atoms or cantilevers, designs a quantum 
interference that eliminates the blue sideband in addition to the carrier transition, thus excluding 
any heating process involving up to one-phonon processes. As a consequence cooling achieves van- 
ishing phonon occupation up to first order in the Lamb-Dicke parameter expansion. Underlying this 
scheme is a combined action of two cooling schemes which makes the proposal very stable under 
fluctuations of the physical parameters such as laser intensity or detuning, making it a viable candi- 
date for experimental implementation. Furthermore, it is considerably faster than existing ground 
state cooling schemes, overcoming one of the limitations of current quantum information processing 
implementations. Its suitability as a cooling scheme for several ions in a trap or for a cloud of 
atoms in a dipole trap is shown. This work extends the description of the "robust cooling scheme" 
introduced in 



I. INTRODUCTION 

Processing quantum information commonly requires 
the initialization of a quantum register to a given pure 
state. In the case of trapped particles such as ions or 
atoms as well as for nanomechanical oscillators, this is 
usually achieved by means of ground state cooling. The 
success of quantum computational tasks, quantum infor- 
mation processing or quantum logic spectroscopy tech- 
niques are directly affected by the accuracy to which 
ground or initial states can be produced. Bose-Einstein 
condensation is reachable only after laser pre-cooling, 
and the very prolific field of quantum simulations relies 
on cooling schemes in many of its implementations. 

The use of laser cooling schemes [2-4] has proven effec- 
tive and useful in most of these fields, with the relevant 
exception of Bose-Einstein condensation, where many ef- 
forts have been directed towards the achievement of an 
all-optical means of obtaining quantum degeneracy [5|. 
Whether the particle is free or bound by an external po- 
tential dictates a fundamental distinction among differ- 
ent treatments. The idea underlying Doppler cooling [6| 
for free particles is related to sideband cooling 0, [sf for 
bound particles and similarly dark state cooling for free 
particles [9] has its counterpart for trapped ions p^of . 

At present sideband cooling is the method of choice 
for ground state cooling of trapped ions. The red side- 
band is preferentially addressed by detuning the laser 
light by the value of the trap frequency v. This requires 
that the linewidth of the transition V allows for the nec- 
essary resolution. Hence this scheme only works in the 
so called strong confinement limit T <^ v. On instances 
where this limit is not satisfied, lower effective linewitdh 
might be engineered by means of laser couplings to other 
levels [11, 12]. Off-resonant heating processes (primarily 
carrier transition excitation, followed by blue sideband 
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heating) limit its performance both in terms of cooling 
rate and final temperature. Its cooling rate is deter- 
mined by the effective linewidth of the optical transition 
Fe// and the coupHng strength of the laser light to the 
electronic levels, corresponding to the Rabi frequency Vt 
times the Lamb-Dicke parameter 77. The minimum reach- 
able phonon number is limited by in the case of very 
low driving. 

The detrimental effect of the heating associated to the 
carrier transition can be overcome by means of dark 
state cooHng schemes. Here, the destructive interfer- 
ence that generates Electromagnetically Induced Trans- 
parency (EIT) [13] cancels the resonant absorption. In a 
three level lambda system, the Raman coupling dresses 
the atomic states giving rise to one dark state and two 
excited states. By adjusting the detuning A and laser 
intensity Vt correctly, an effective coupling of the dark 
and the longest lived excited state can be achieved with 
a detuning equal to the trap frequency [13, • The in- 
volvement of the dark state ensures the cancellation of 
the carrier transition, and the detuning adjustment en- 
hances the red sideband transition with respect to blue 
sideband excitations. All in all, final ocupation numbers 
proportional to can be achieved, while the rate 

scales as rf^^ thus beating sideband cooling for large 
detunings. 

In the case of the Stark shift (SSh) cooling method 
[l6|, it is assumed that the Raman coupling present in 
the EIT scheme does not affect the mechanical degrees 
of freedom. Instead, this role is taken over by a direct 
and resonant coupling of the two lowest lying states. The 
dark and bright states with respect to the Raman cou- 
pling are precisely the dressed states of this new cou- 
pling, and its Stark shift can be adjusted through the 
coupling intensity [17]. If adjusted to the value of the 
trap frequency, a red sideband coupling is favored. This 
effectively transfers the mechanical energy to the bright 
state and then it is dissipated through its coupling to 
the excited state. Higher values of laser intensities can 
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be applied with this proposal, allowing for an effectively 
faster operation of the scheme. 

The limiting factor on both EIT and SSh coolings is 
the heating associated to the off-resonant blue sideband 
which, after the carrier transition, is the only heating pro- 
cess left up to first order phonon processes. For large F it 
is not possible to neglect the effect of the blue sideband, 
and cooling efficiency might be affected. 

In [1] it has been shown that the blue sideband tran- 
sition can be effectively canceled by the combination of 
both EIT and SSh cooling schemes. Both schemes are 
using the same dark state but different mechanisms, so 
that their combined effect on the same system is that of 
effectively canceling the blue sideband in addition to the 
carrier transition. Furthermore, the condition for cooling 
depends on the ratio of the Lamb-Dicke parameters of the 
couplings involved (rather than on the Rabi frequencies) 
thus making it a more experiment-friendly scheme. 

In this paper we set out in detail the mechanism under- 
lying the proposed scheme, with particular emphasis on 
possible systems for experimental implementation. The 
paper is organized as follows. In section [IT] an explana- 
tion of the mechanism of the cooling scheme is presented, 
while the detailed explanation of the theoretical treat- 
ment is introduced in section UTTl An analysis of the prop- 
erties of the scheme begins with section IIV| where the 
robustness with respect to fluctuating experimental pa- 
rameters is explained. The optimal cooling rate is found 
in section |Vl In section IVll two possible experimental im- 
plementations are set out, and the possibility of cooling 
the motion on more than one axis with the same scheme 
is presented in section IVlIl Section IVlIII deals with the ef- 
fect that phase mismatch of the cooling fields would have 
on the efficiency of the scheme. The last two sections are 
devoted to the study of the appHcability of this scheme 
to many body systems. On the one hand, in section HXl 
the situation is considered where several ions sit on the 
same trap. On the other hand, in section |X] a cloud of 
dipole trapped atoms is assessed. 



II. MECHANISM 

Our proposal is designed for a 3-level lambda system of 
mass m which is trapped in a harmonic well of frequency 
u. This can be an accurate model for an ion in an elec- 
tromagnetic trap or an atom in a deep dipole trap. The 
three levels are coupled by means of an electric dipole 
interaction with running waves as shown in figure [T] and 
expressed in the following Hamiltonian: 

H = ub^b + uje \e) (e| + uj^ |t) (t| + oJi \i) {i\ 

+ l^AC^i^'^^ COS {ujAt + kAx) ^ 
+ rtAcr^''''^^ COS {uAt - kAX) 
+ VLbct^x'^^ cos {ujBt + ksx) , 




FIG. 1. A 3-level electronic system consisting of |t) and |^) 
and an excited state |e) which dissipates energy at rate F. The 
lower levels are coupled to |e) by a pair of Raman beams under 
a detuning A whose Rabi frequencies are and wavevector 
projection on the cooling axis of opposite sign (this may 
be achieved, e.g., by counterpropagating beams). |t) and |^) 
are directly coupled by a beam of Rabi frequency Q.b and 
wavevector projection ks- 



where cfx' = |<^) (6| + |6) (a|, VLa.b are the Rabi frequen- 
cies of the respective laser couplings, C(Je,t,i are the energy 
of the respective levels and oOa ^ are the frequencies of 
each laser. The metastable state |t) and the ground state 
1^) are considered to be infinitely lived, so the only dis- 
sipative state is |e). 

This time dependent Hamiltonian can be simplified by 
moving to an interaction picture with respect to the ap- 
propriately detuned energy terms of the internal degrees 
of freedom. This requires the definition A = cjg — ~ 
ijOa = oje — oj^ — and setting uj^ — uj^ = ujb- High fre- 
quency terms are dropped under a rotating wave approxi- 
mation. It is useful to express each wavevector projection 
in terms of its corresponding Lamb-Dicke parameter fol- 
lowing the definition rj = ^/cxq, where xq = is 
the zero point motion of the oscillator. Throughout the 
paper we consider the trapped system to be in the Lamb- 
Dicke regime, in which 77a,^s ^ 1- In physical terms, 
this implies that the recoil energy gained in each photon 
emission is much smaller than the necessary energy to 
move one level in the trap, therefore processes involving 
phonon creation or annihilation are realized with small 
probability. An expansion of the Hamiltonian up to sec- 
ond order phonon processes is hence justified. Expressing 
the subspace spanning , with the basis formed by 
l+> = 71 (It) + and I-) = ^^(It) - i;)) the Hamil- 
tonian can be split into the following 4 terms: 



H = Htr+Hi„t + VEIT + VsSh, (2) 
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with 



A. Derivation of the Steady State 



Htr 


= vh%, 




= A|e) (e|, 


Veit 








VsSh 






V2 





(3) 



where ai^'^^ = \a) {a\ - \b) 4""'^^ = ^ !«) (^1 - ^ 1^) («l 
and ^ = -^{b -\- b^). Veit and Ks^/i are also refered to 
as the A and B couphngs respectively. This Hamiltonian 
can be regarded as a combination of those describing the 
reversible dynamics of EIT and SSh coolings as expressed 
by 



Heit = ^b^b + A |e) (e| + Q.a 



,(e,+) 
>x 



(4) 



and 



Hssh = + A|e) (e| 



Ha 
7=2 



.(e, + ) 



Ob 



(5) 



In the EIT case, the cooling transition is |— ) ^ |e). 
There is no direct carrier transition driving it, only a 
combination of a blue and a red sideband expressed by 
the last term in (|4]). After diagonalization of the second 
and third terms, the values of VLa and A can be adjusted 
so that the detuning to one of the resulting dressed states 
matches the value of the trap frequency In this way, the 
only process under resonance will be the red sideband to 
the resonant dressed state, and all other opto-mechanical 
processes will be off-resonant, including the two photon 
processes coupling W) and |^). In the limit of large A, 
though, these two photon processes would gain in rele- 
vance. 

The situation is similar in the SSh case. This time 
the cooling transition is |— ) ^ |+). The carrier coupling 
between |t) and \\) transforms into a Stark shift in the 
1+) and |— ) picture. This Stark shift can be adjusted 
to the value of the trap frequency so that again the red 
sideband is in resonance and the blue sideband is not. 

This presents us with a Hamiltonian capable of per- 
forming two cooling schemes that apply to the same sys- 
tem and that operate by means of independent mech- 
anisms. Both can be combined so that the only heat- 
ing process remaining in both cooling transitions, namely 
the blue sideband, is canceled. In a master equation de- 
scription spanning only up to one phonon processes, this 
amounts to only cooling processes being present; there- 
fore the steady state of the system must be one with zero 
temperature, i. e. a pure state. 



It is possible to derive the form of the pure steady state 
of the system by analyzing the splitting of the Hamilto- 
nian into EIT and SSh parts. What we want to find out 
is if there exists a steady state of the system that is pure 
and under which conditions this is so. 

A steady state is such that the master equation van- 
ishes. The coherent contribution of the Liouvillian van- 
ishes for any eigenstate of the Hamiltonian, while the 
only atomic state involving incoherent dynamics is the 
excited state. This shows that a plausible instance of a 
pure steady state would be an eigenstate of the Hamilto- 
nian with vanishing overlap with the excited state. 

The only part of the Hamiltonian coupling the state 
|e) with the rest of the Hilbert space is the third term 
in eq. [3] {Veit) related to the EIT part of the cooling. 
Therefore, it is a condition that the eigenstate be dark to 
this term, i.e., Veit \^) ss ~ ^- Since the system is within 
the Lamb-Dicke regime, perturbation theory is sufficient 
to solve it. Notationally, the Lamb-Dicke parameter or- 
der is indicated by a bracketed superscript. At leading 
order the equation becomes cr^'^^ \^)^ss ~ ^' which im- 
plies = | — ) 1^)^, with |(/))^ an undetermined state 
for the mechanical degrees of freedom. For the first order, 
crt^^ + V2VA(Ti'''~^q\'^)fJ = 0. This determines 
that \^)i'J = -zV2vA 1+) where \^^)^ = q \^)^. 

The eigenvalue equation for the rest of the Hamilto- 
nian can now be applied to fully determine the state. In 
the zeroth order, the only term involving the mechanical 
degrees of freedom is the frequency term Htr = ^b^b. As 
a consequence, the mechanical part of the zeroth order 
eigenvalue must be a Fock state, so l^)^^^ = |— ) |^) and 
a/2 |(/>0^ = \/ri |n — 1) + -\- 1 |n + 1). The next order 
equation determines the value of n. Since V^^^ l^)is^ is 
proportional to l^)^^^ the eigenvalue equation requires 
that {Htr + Vs^i) also be proportional to 1^)1^^ 

This will happen for |^0m ^ Fock state as well. The only 
value of n for which both |n) and q \n) are Fock states is 
n = 0. 

The state l^)^^ has been fully determined as |— ) |0) — 
^Va 1+) Whether this is indeed an eigenstate of the 
Hamiltonian can be guaranteed perturbatively, but the 
eigenvalues at each perturbative order have to be equiv- 
alent. Since the eigenvalue obtained at zeroth order is 
and the one at first order is + ^(1 - ^), the 
following condition arises: 

r]B _ ^ 
r]A 



2. 



(6) 



This ensures that l^)^^ is the unique eigenstate of the 
Hamiltonian that is a steady state of the master equation 
at the same time. This condition balances the values of 
the Lamb-Dicke parameters with respect to the value of 
the laser intensity and is shown below to be very robust 
to fiuctuations. 
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B. Double Interference 

This scheme couples the double transition | — ) ^ |e) 
and |— ) ^ 1+) to the mechanics of the particle to gener- 
ate cooling. It has been shown that the state |— ) is de- 
coupled from the rest of the Hilbert space at zeroth order 
in the Lamb-Dicke expansion. This constitutes the first 
interference of the system, canceling any carrier heating 
in the cooling transition and isolating | — ) as the zeroth 
order electronic steady state. The electronic degrees of 
freedom unperturbed by the motion are described by the 
Hamiltonian 



int 



(0) 



EIT 



(0) 



SSh 



i=e,+ 



(7) 

where, taking the level | — ) as the origin of energies, Se 
corresponds to A + ^ and ^+ is just Qb- The Qa cou- 
pling dresses the excited and bright states into \Di) and 
\D2) in the form 



1+) =cos6>|L>i) + sini9|L)2), 
|e) =sinl9|Di) -cos(9|D2), 



where 



tan^ : 




(8) 
(9) 



(10) 



The eigenfrequencies of the states are Sjj^ ^ = 

S^^S^T \J'^^\ + (^^e - ^+f^ /2. Both dressed states 

are dissipative due to their finite overlap with the excited 
state. Their energy uncertainty is Lorentzian-shaped but 
each has a different width proportional to the overlap to 
|e). In terms of the unperturbed states, though, the cou- 
pling to the continuum is the addition or subtraction of 
two Lorentzians. The case of the excited state is particu- 
larly interesting, since the change in sign makes it couple 
with positive amplitude for high energies and with a neg- 
ative one for small ones. This together with the always 
positive coupling of the bright state can be combined to 
generate an energy-selective interference. An excitation 
from the dark to the excited state will interfere with that 
to the bright state for a particular value of energy. This 
value can be selected to be precisely that canceling the 
blue sideband by correctly adjusting the experimental 
parameters, as is derived below. 

Because the excited state couples to a continuum of en- 
ergies, a complete description should involve the modes 
of electromagnetic radiation to which it couples, so that a 
state of the system is described by the infinitely degener- 
ate energy subspaces {|z,n)} with i any of the electronic 
levels and n the number of photons in any of the infinitely 
many modes of the environment, from where degeneracy 
arises. This notation is not to be confused with that 
of previous section, where the second digit corresponded 
to the phononic number. Assuming the radiation modes 



other than those of the lasers are empty, the Hamiltonian 
above would involve the states |z, 0) (which are not degen- 
erate and can therefore be redefined An additional 
two couplings |e) ^ |+, 1) and |e) ^ | — , 1) complete the 
description of the system. 

The manifold {|+) , |e) , {|+, 1)}, {|-, 1)}} together 
with the corresponding couplings can be diagonalized, so 
that this would define a new continuum \k) with eigen- 
value k. It is then possible to derive the relationship (see 
Appendix 



(e|fc) = ^(fc-f^B)(+|fc) 



(11) 



This introduces a change of sign of the overlap for en- 
ergies below that of the state |+,0), and the crossing is 
known as a "Fano resonance" . This change of sign is es- 
sential for the emergence of the interference in the blue 
sideband interaction. 

Indeed, the blue and red sideband interactions can be 
expressed as 



(1) 



EIT 



(1) 



SSh 



(12) 

To zeroth order we assume the electronic state of the 
system to be located always in | — ). The interaction can 
hence be simplified to 



(13) 



Note that the notation including the photons has been 
dropped. The levels |e) and |+) have been shown to 
be a superposition of free energy states \k). Energy 
conservation will force the system to preferentially per- 
form the transitions \k = v){—\a (the red sideband) 
and \ — h') {—\a^ (the blue sideband), as are depicted in 
Fig. ([2]). The strength of these two transitions for the 
current interaction operator are 



(14) 



and, similarly, for the blue sideband, 

(-H+)|-^)(-|a. (15) 



-r]A {^B + ^) 



The blue sideband can be made to vanish if the equation 
Va {^B -\- ^) — = is satisfied. This condition 

coincides with Eq.®, which guarantees that the steady 
state is a pure state. This indicates the equivalence of 
both requirements. 



C. Fano Resonance 

From equation (pTj) it is clear that there exists an en- 
ergy eigenstate from the continuum \k) that is orthogo- 
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tally accessible, making this option immune to the effects 
discussed in section IVIIIi 
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FIG. 2. Representation of the two first order processes in the 
system. In red (solid line) is the red sideband, which decreases 
the mechanical energy by one level while the internal degrees 
of freedom get excited to an energy of The opposite process, 
in blue (dashed line), doesn't occur because the probability 
amplitude of the path through |e) and through |e) interfere 
with each other. 



III. ANALYTICAL TREATMENT 

Following the procedure in [is"-^ the system under 
study is described using a master equation formalism, 
so that both the coherent dynamics and the dissipative 
nature of the excited level can be accounted for: 



dp 
~dt 



(16) 



where p is the state of the system involving both the 
internal and the external degrees of freedom. The su- 
peroperator is a Lindbladian for the two dissipative 
channels: 



where aj^k = \j) (^|, the same rate 
assumed for both channels and 



7e,t 



F is 



nal to |e). This is the phenomenon known as Fano res- 
onance, and it offers an alternative way to implement a 
double interference in the system. For this, the coupling 
of the states |t) and ||) needs to be zero at first order 
in the Lamb Dicke expansion (i. e., tjb = 0) and the 
Rabi frequency needs to satisfy the resonance condition 
Qb = — ^- The nature of the interferences is therefore 
different as in the general case. For the Fano Resonance, 
the interference arises between the coupling of |e) and 
the continuum on the one hand and the coupling of |+) 
and |e) on the other hand. In the general case, the in- 
terference is between the laser couplings |— ) ^ |+) and 
|— ) ^ |e). Although the principle behind both alterna- 
tives is different, condition ([6]) contemplates both, since 
it gives the correct condition for Qb when tjb = 0. 

The relevance of this alternative is paramount, since 
this might simplify considerably some experimental im- 
plementations. One of the challenges associated with the 
implementation of the scheme presented in this paper 
might be coupling |t) and |^) with a Lamb-Dicke pa- 
rameter that is on the order of an optical one, since the 
energy difference between them is not of optical range. In 
section IVll some options to generate optical Lamb-Dicke 
parameters are presented, involving two photon processes 
or magnetic gradients. The B coupling in the Fano reso- 
nance alternative doesn't need to involve the mechanical 
degrees of freedom, i.e., its Lamb-Dicke parameter can 
be negligible. Hence, non-optical means of coupling can 
be exploited, like microwave radiation or Zeeman effects. 
This makes the scheme more accessible where geometrical 
adjustments of the Lamb-Dicke parameters (described in 
section IVip are not possible. In addition, phase locking 
a microwave coupling to the laser coupling is experimen- 
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dsV^(5)e^^^'^^>e-^^^'^^" (18) 



accounts for the momentum transfer of Hke^i in the event 
of a photon emission concomitant to the electronic decay 
from level |e) to level W{s) = |(1 + s^) is the an- 
gular distribution for a spontaneous emission of a dipole 
transition. 

In the Lamb-Dicke approximation the Lindbladian can 
be expanded like the Hamiltonian, so that the first three 
terms of the whole Liouvillian C would appear 



A|e)(e| + ^cx(-+) + ^cx(+'-), 
^^oe{p) +>Co/(p), 



VBCT^y^' ^ ]q,P 



iLi^2 (e,+) , (+,-) 



f^P 



^ Je^iVfcTi^ei'^qpq ^ f P ^ pf)cre,i 



(19) 

where rji is the Lamb-Dicke parameter corresponding to 
/ce,i. The zeroth order Liouvillian has been split in a part 
acting only on the external degrees of freedom Cqe and 
another one acting on the internal degrees of freedom 
jCqi to emphasize a lacking coupling among both at this 
order. 
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In the attempt to approximate a solution to the mas- 
ter equation, the theoretical approach in [11] is followed. 
Since the internal and external degrees of freedom are 
not coupled in the zeroth order Liouvillian, any steady 
state is separable, so that it can be expressed as a tensor 
product of electronic and mechanical states. Following 
the analytical procedures applied in the previous section, 
any zeroth order steady state is of the form |— ) \n) with 
n = 0, 1, .... The leading order Liouvillian has hence 
an infinitely degenerate subspace of steady states, i.e., 
of eigenvectors with an eigenvalue equal to zero. This 
eigenspace is connected by Ci and £2 to the rest of sub- 
spaces. To the extent that the cooling rate is slower 
than the dynamics of internal degrees of freedom and 
the trap frequency, it is possible to disregard the rest of 
the Hilbert space and focus on an effective master equa- 
tion describing the dynamics of the phonon populations. 
This is achieved by projection of the global master EqIT6l 
to the null-eigenvalue subspace up to the second order. 

^ = [VC2P + P£i(-£o-^)£iP]p, (20) 

This can be used to derive a rate equation for the popu- 
lations only: 

±{n) = -{A.-A+){n)+A+, (21) 

where and A_ are the heating and cooling rates re- 
spectively and they are a function of the system parame- 
ters with the condition A-^(iy) = A-{—u). This offers an 
analytical prediction for the steady state mean occupa- 
tion number (n) and for the rate of the cooling process 
W within the range of validity of the approximation: 

/ \ - 

(22) 

W = A_ - A+. 

For the present scheme, the derivation of the rates gives 
the result 

= [2r]A{jy + Qb) - VB^B^/V, (23) 

with n\rV = 4r^{iy^5^)^^4{Q\/2-{iy^S^){u^6e)f. 
This rate will vanish exactly under the condition ex- 
pressed in Eqini confirming the reasoning of Hamilto- 
nian and coupling interference presented in section HIl 
Vanishing heating rate indicates that the only remaining 
heating process in the system, namely the blue sideband, 
is being effectively canceled. This also implies that the 
cooling rate W is for that case equal to A- . 

IV. ROBUSTNESS 

The destructive interference between EIT and Stark- 
shift contribution is crucial for understanding the robust- 
ness of the scheme under fiuctuating parameters. If the 



0.00 




FIG. 3. Expected value of occupation (n) as a function of the 
Lamb-Dicke parameter tja and of variations around the opti- 
mal Rabi frequency Q^, where it is shown how fluctuations 
have a worse effect as r/A vanishes, that is, as the Stark-shift 
cooling regime is approached (77s = 0.4, Qa = O.lu). 

Rabi- frequencies deviate from Eql6]by Aft^/B^ the final 
population is affected by 

(n) oc (An^)4(Al^B)2, (24) 

in contrast to the second order dependence as is usually 
the case. In other techniques such as SSh or EIT the res- 
onance conditions do not necessarily match a population 
minimum and fluctuations can be very signiflcant. As is 
exemplifled in Figl3l operation of Stark-shift only or EIT 
only cooling schemes suffer more from Rabi frequency 
fluctuations than the combination of both in the "robust 
cooling scheme" . This can act as an experimental protec- 
tion to the performance of the cooling scheme, providing 
more certainty to the theoretical predictions also in a 
more realistic frame of fluctuating laser intensities. 

Indeed, actual experimental application of EIT cool- 
ing yielded flnal population values of order 10~^. The 
theoretically predicted value is lower of about an order 
of magnitude, indicating that a realistic prediction of the 
flnal result should take into account the uncertainty in 
the experimental parameters. Our analytical model pre- 
dicts flnal populations of the order of 10~^ even for Rabi 
frequency fluctuations of about 2%, leaving thus ample 
room for experimental improvement. 

The robustness of the system can be improved signif- 
icantly if the Rabi frequency of the B coupling is set 
as close as possible to the value of the trapping fre- 
quency and Eql6] is accordingly met. In this circum- 
stance the contribution to fluctuations from the term 
{AQa)'^{AQb)'^ almost vanishes and the next term in 
the Taylor expansion takes over, thus yielding 

(n) oc (A^^A)^(A^^B)^ (25) 

which ensures excellent stability and hence makes it ar- 
guably advisable, although not mandatory, to set the 
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scheme to meet VLb = ^- As is shown below, this value 
is not optimal for cooling rate considerations, so a com- 
promise between both is necessary. 



0.07 r 



V. COOLING RATE 

The analytical prediction for the cooling rate derived 
above reads 



W 



(26) 



This rate corresponds to the absorption probability at the 
trap frequency since, under condition (j6j), the emis- 
sion at that frequency vanishes. The absorption 
spectrum has a peak for each of the dressed states \Di) 
and 11^2), so the rate will be highest when the broader 
of the two resonances matches the trap frequency. The 
conditions = ^ corresponds to 



(27) 



and the inequality 5+ > 6e > so that the resulting rate 
expression is 



W 



2 • 



(28) 



The vahdity of this expression is constrained to the per- 
turbative treatment applied here. The internal dynamics 
governed by ftifji must always be slower than the inter- 
nal dynamics. Hence, it breaks down for cases where 
ft A ^ ^e,^+, since then it is \Di) the broader dressed 



state. The singularity S-^ = u is resolved as 



al- 



though it corresponds to vanishing and is then of 
no interest. The validity of the expression is explored 
in Fig.(j4]) for a particular set of parameters, where it is 
shown how the analytical calculation matches the numer- 
ical simulation only for small values of Qa- 

Placing the resonance at the value of the trap fre- 
quency imposes one condition on the parameters, but 
there are still enough degrees of freedom to lock the po- 
sition of the second peak. The first order laser-motion 
interaction Eq.([T2j) generates two phonon processes when 
second order perturbation theory is invoked. In order to 
enhance this effect, one can place the peak at two times 
the trap frequency. The conditions 6^2 = ^ and S^i = 2v 
can then be simplified to 



(29) 
(30) 



This and the fact that Aj^ vanishes in this scheme allows 
for faster cooling rates as compared to schemes where 
this is not the case. As an example, the plot in Fig.(j5j) 
compares the cooling rate as a function of VLa for our 
scheme and EIT cooling. Because SSh is only slightly 
slower than our scheme makes its curve overlap and it 




FIG. 4. Comparison of the cooling rate W obtained by means 
of numerical (solid blue line) or analytical (dashed red line) 
calculations as a function of the rabi frequency (F = 15z/, 
Q.B — 1.3z^, riB — 0.1). For large values of the analytical 
result fails to predict the numerical values. 




FIG. 5. Comparison of the numerical cooling rate W for our 
scheme (solid blue line) and EIT cooling (dashed red line) as a 
function of the rabi frequency Q^a (F = 15z^, Q^b — 1.3z^, rfB — 
0.1). The current scheme has in general a better performance 
already for smaller Rabi frequency values. 



can't be shown in the plot. Optimization of the values of 
both Rabi frequencies for F = 15z^ and = 0.1 can pro- 
vide cooling rates of up to 0.06z/, whereas higher values 
for the Lamb-Dicke parameter makes it possible to accel- 
erate the scheme to just one order of magnitude below 
the trap frequency. 
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VI. EXPERIMENTAL IMPLEMENTATION 

This proposal can be experimentally implemented in 
a number of different ways. Although the scheme ap- 
plies to very general three level systems, we assume the 
ground level transition to be of microwave or rf order and 
the excited state to be at an optical distance above the 
ground level. The A coupling can hence be performed 
with two optical Raman laser beams, and the question 
remains then of how to perform the microwave coupling 
B. The scheme condition Eql6] lays a constrain for the 
Lamb-Dicke parameter of coupling B, which has to be at 
least double that of coupling A. Another option is avoid- 
ing mechanical coupling of the ground state interaction 
{rjB = 0), following the "Fano resonance" alternative of 
the scheme. This simplifies the scheme in that a mi- 
crowave coupling would be enough and locking it to the 
laser is experimentally accessible, making it resistant to 
phase effects discussed in section IVIIII Two viable im- 
plementations are taken under consideration below, al- 
though they don't exhaust the range of possibilities. 



A. Magnetic gradients 

The direct driving of a microwave transition has an as- 
sociated Lamb-Dicke parameters which is orders of mag- 
nitude smaller than that of an optical frequency. In order 
to increase its val ue, a magnetic gradient can be applied, 
as is set out in [U, [22|. In this system the magnetic 
gradients create a coupling of the following type: 

Aa,(6 + 6t), (31) 

where A is proportional to the gradient of the magnetic 
field and the two level system is driven using a microwave 
of the form QdO'x cosudt, where ftd corresponds to the 
Rabi frequency and Ud to the angular frequency of the 
driving wave. After a Schrieffer- Wolff transformation the 
resulting Hamiltonian is exactly as in Eql3]when the Rabi 
frequency is replaced by Qd and the Lamb-Dicke param- 
eter is replaced by ^ . The range of practically achievable 
values is comparable to that of optical transitions [21]. 

This scheme can be especially useful to cool nano scale 
resonators, by using the setup described in [23]. In this 
setup an NV center is coupled to a diamond cantilever, 
the coupling is performed by magnetic gradients result- 
ing in the same Hamiltonian as described above. In can- 
tilevers the speed of cooling is very important due to 
the finite Q value, which is a central factor limiting the 
attainable final temperatures at present. The high cool- 
ing rate achieved by the described scheme will result in 
lower final temperatures bringing us closer to the goal of 
reaching the quantum regime in cantilever systems. 




FIG. 6. Physical realization of and Q.b couplings, for 
cooling in the trap axis, as two Raman pair of beams char- 
acterized by their respective Rabi frequency Q^a and and 
their Lamb-Dicke parameters rj^ and rjp. 



B. Highly detuned Raman beams 

Alternatively, the ground states can be coupled by op- 
tical means, thus directly ensuring a Lamb-Dicke param- 
eter of similar order as for coupling A. In order to do 
so, Raman beams with large single-photon detuning A' 
are envisaged to couple levels | t) and | ^). As is shown 
in Figini they have independent Rabi frequency Vtp and 
Lamb-Dicke parameter rjp. The large detuning effectively 
decouples it from the excited state, unlike coupling A. 
Hence, adiabatic elimination of the upper level is possi- 
ble and it yields the relationships between our effective 
parameters VLb and r]B and the physical values and 
Tjp. Its derivation is presented in appendix [Bl and it is 
found that VLb = ^p/A' and tjb = '^rjp for sufficiently 
large detunings (A' <C v) and in the Lamb-Dicke regime. 
Meanwhile, the beams for the A couplings are tilted an 
angle 6 with respect to the trap/cooling axis, so that 
Va = r]'j^ cos 0. This provides us with an additional geo- 
metrical degree of freedom to adjust the ratio r]B to r]A- 

Even though the detunings of each of the couplings 
are different, the fact that we are dealing with optical 
frequencies justifies 77^ — ^p, and under this assumption 
— = ^-n. By relating this to the resonance Rabi fre- 

7] A COS U ^ " 

quency, cos = 77^^^ , it is shown that the available angle 
range ^ = [0, |^] will span all possible values of Qb- If we 
choose the optimal point for the fiuctuations about the 
value VLb = the condition becomes tjb/va =4, which 
corresponds to a layout where = 60^. A compromise 
with the cooling rate has to be reached, though, so that 
the optimal angle will be close but not exactly 60^. 

It may happen that the trap has very limited optical 
access so that is constrained to other non-optimal val- 
ues. There is a way around this on the grounds that 
A' still constitutes a useful degree of freedom. Since it 
has no upper constrain, the B coupling can be performed 
with frequency doubled beams so that r]p 2r]A and the 
optimal constrain is achieved already at ^ = 0. 

Additionally, in cases where neither the geometrical 
configuration nor the detuning can be optimized, the 
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aforementioned robustness of the scheme ensures excel- 
lent performance. Whether the cooling rate or the final 
temperature need to be improved, appropriate set of pa- 
rameters can be found by numerical means. Usual exper- 
imental layouts of ion traps use vacuum chambers with 
windows at 22.5^ and/or 45^ from the trap axis, which 
generally allow for an angle range of about ±10^. As- 
suming A' can't take values on the optical range, and 
taking 45^ as an operating value, the Lamb-Dicke pa- 
rameter quotient becomes tjb/va = 2^/2. Eql6]can still 
be observed by adjusting Qb^ and Qa = 0.6i^ and A 
optimizes the cooling rate for this condition. This final 
result can be improved depending on the particular val- 
ues of the transition linewidth F. On the contrary, if the 
cooling rate is to be enhanced, Eql6]will not be satisfied. 
In particular, for Qa — 0.4i^, — OAbu with 5% fiuctu- 
ation and A —2v the population can still be as low as 
10~^ while having a cooling rate faster than that of EIT 
cooling. Taking into account the fact that angles up to 
55^ are accessible the cooling rate can still be improved 
by up to two orders of magnitude. 

Actually, the implementation picture of a double pair 
of Raman beams makes both couplings totally inter- 
changeable, allowing for further physical intuition. While 
a configuration of the two A beams perpendicular to the 
trap axis represents the Stark-shift cooling (disregard- 
ing Eql6] and tuning VLb = a value of ^ = and 
the corresponding Rabi frequencies tuning is equivalent 
to EIT cooling. Any intermediate configuration realises 
an instance of Robust cooling. This picture is also very 
intuitive to understand how fiexible the cooling scheme 
is regarding the cooling of the radial modes, since the 
addition of only two more laser beams in a suitable in- 
clination with respect to the trap axis can cool down all 
three dimensions in the same Robust regime. 



VII. MULTIAXIAL COOLING 

So far one-dimensional cooling has been discussed. 
Nevertheless, the fiexibility of the implementation of fig. 
(|6]) facilitates addressing more than one direction in the 
trap, so that multiaxial cooling can be achieved without 
adding any other beams to the scheme. In particular, 
two independent azimuthal angles constitute the neces- 
sary degrees of freedom to adapt the cooling for two per- 
pendicular directions. The question remains whether the 
third direction will as well verify condition (|6|) and to 
which extent. This is in general not going to be the case, 
but even in this situation detuning degrees of freedom 
can be used to implement Raman sideband cooling to 
that axis. 

In the case where traps have very limited access it is 
always possible to find an axis where cooling takes place. 
When two independent dimensions share the same trap- 
ping frequency it is always possible to find an axis for 
which the cooling condition is fulfilled. This is usually 
the case for the both radial modes in an ion trap. As 



Radial plane 




,,-Oprtirnal Cooling Axis 



FIG. 7. Geometrical definition of angle 0' and representation 
of the projection of the two Lamb-Dicke parameters to the 
optimal cooling axis when 77^ — Tjp. 



long as the projection on the radial plane of couplings 
A and B corresponds to two pairs of counterpropagating 
beams 6>-rotated from each other, it is possible to define 
an axis where the projection of the couplings fulfills the 
cooling condition. This axis can analytically be charac- 
terized by the definition of a new polar angle 0' as shown 
in Figlll 

The ratio of Lamb-Dicke parameters A and B on any 
axis 9' of the plane is 



r]B 2co^6' 
^ " cos((9 + 6') 



(32) 



The newly defined axis will fulfill the optimal cooling 
condition (|6]) for the angle 



6' = arctan 



(l^B + ^) cos6> — VLb 
(ytB + ^) sin6> 



(33) 



For any angle in the range (0, ^) there always exists 
a solution in the range (0, Furthermore, it can be 
shown that the sum of both 0' and is never larger than 
7r/2, so that geometrical considerations can be restricted 
to one quadrant. 

Only geometrical degrees of freedom have been ex- 
hausted in this description. For traps with limited ac- 
cess also detuning degrees of freedom can be exploited in 
order to guarantee cooling in all three spatial directions. 



VIII. PHASE EFFECTS 

In view of an experimental implementation it is impor- 
tant to determine whether beam phase mismatches can 
be relevant to the performance of the scheme. Surpris- 
ingly enough, out of 4 beams used in the implementation 
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FIG. 8. All possible phase mismatches expressed with respect 
to the A coupling (red) between |t) and |e). 



proposal a single phase mismatch can account for all the 
effects on the performance. To show this, it is helpful 
to define all possible phases, as in figlH with respect to 
the beam of coupling A between |t) and |e). The sys- 
tem would be oblivious to phase mismatch in the case 
where /3 = and a = (j), i.e., to a global phase difference 
between beams coupling different ground states. This 
is so because one can always redefine the corresponding 
ground state so that it absorbs the same amount of phase 
in both beams. 

Out of the pair of Raman beams performing the B cou- 
pling, only the difference of its phases survives in the ef- 
fective microwave coupling between ground states. This 
effect leaves room for a redefinition of their phases so 
that a ^ {a — (3) and (3 ^ is rather chosen. After 
the redefinition, the phase difference between both legs 
remains the same, i.e., the effective Hamiltonian is un- 
changed while the initial picture is simplified to a 2 phase 
system. 

Now only the two beams coupling |^) to |e) have a non- 
zero phase. The redefinition ||) g-<a-/5) ||^^ leaves 
the beams for the B coupling phaseless while the left A 
beam is the only one carrying a phase. As a conclu- 
sion, the phase originally defined as (j) accounts for all 
the effects that laser beam dephasing could have on the 
system. 

Numerical simulations of the effect of a given phase 
for the final achievable temperature and the cooling 
rate are shown in fig. [9] and [TOl The truncation values 
puts a limit to the accuracy of the results when they 
get closer to the maximum occupation number. Thus, 
for the final temperature results, any value of the final 
occupation over 5 is to be taken as unreliable. If the 
incidence angle of the A coupling is the same as for the B 
coupling it is very likely that both beams follow a much 
more similar path and that phase effects can be controlled 
much better. 

In addition, for the "Fano resonance" alternative of im- 
plementation, these effects are going to be less relevant. 
It is possible to phase lock a microwave and a laser field 
so that phase mismatches are reduced to a minimum. 



0.12 
0.1 
0.08 
cT 0.06 
0.04 
0.02 




-8.1 -0.05 0.05 0.1 



(f) (rad) 



FIG. 9. Effect on final temperature of dephasing of one of the 
beams (Qb — 5z^, — 400i/, F — 15i/, r/B — 0.1). 
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FIG. 10. Effect on final temperature of dephasing of one of 
the beams (Qs 5z/, Qa 400z/, F 15z/, r]B 0.1). 



IX. MULTI-MODE COOLING 

The cooling scheme has also been tested for an ion 
chain using Monte Carlo simulation [25|. The ro- 
bustness of the scheme implies a wide range of oper- 
ational Rabi frequencies or, in a different perspective, 
a wide range of trap frequencies for a given Rabi fre- 
quency. Thus, a particular central mode frequency can 
be addressed so that also the neighboring modes benefit 
from the cooling. Numerical test have been performed in 
a multi-mode environment with up to 3 ions and promis- 
ing results have been obtained. 

Quantum jumps method, also called Monte Carlo 
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FIG. 11. 3 mode cooling where mode 2 is addressed. While 
the cooling is fast and to a low temperature in all modes, 
neighbouring modes 1 and 3 see a slight decrease in cooling 
efficiency {Qb = z^, = 2.3z/, F = 15i/, tjb = 0.1, with u the 
eigenfrequency of the second mode). 



wave-function method, is a useful tool to simulate the 
performance of our Robust cooling scheme in a multi- 
mode environment. A benefit from this method is found 
mainly for a system with a large number N of states, 
since it is equivalent to a master equation approach to 
an open system while still using the pure state formal- 
ism, thus reducing the number of tracked variables from 
N'^ to N. It is based on evolution with a non-Hermitian 
Hamiltonian and random quantum jumps, followed by 
wave- function renormalization to account for the non uni- 
tary evolution of the non-Hermiticity. The time step St 
that can be used is lower bounded by the correlation 
time of the reservoir and upper bounded by the limita- 
tion r]i5t <C 1, with rji the eigenvalues of the Hamiltonian. 
At the time of simulating, a time step that allows for a 
quantum jump with a probability of roughly 10% is a 
good benchmark. 

A simulation code has been developed that applies this 
method to a system composed by a chain of three level 
ions in a harmonic potential. It has been tested and 
benchmarked with side-band cooling of a single ion and 
then simulations were performed for three level systems 
as instances of EIT cooling, Stark-shift cooling and the 
proposal on this report. Observations led to the con- 
clusion that our proposal is able to transfer the cooling 
ability to neighboring modes much more efficiently than 
other schemes, although as expected cooling rates and fi- 
nal temperatures are not as good as for the central mode. 
An instance of this performance is shown in FigJTTl 



X. ATOM CLOUD COOLING 

If Bose-Einstein condensation could be obtained by an 
all-optical means instead of evaporative cooling, higher 
densities could be obtained, since it would be possible 
to keep most of the initial atomic sample in the trap. 
The longer life of the sample would allow for recooling 
at the end of the experiment since it is a non-destructive 
procedure. 

When, instead of a single atom, it is a collection of 
atoms that is to be cooled, the photons emitted by one 
atom can be reabsorbed by the next one, giving rise to 
heating mechanisms that need to be accounted for in the 
theoretical description. Independently to photon reab- 
sorption, the emission lineshape of an atom is also af- 
fected by the presence of other atoms in what is called the 
phenomenon of superradiance. These processes, known 
globally as dipole-dipole interactions, become significant 
when the distance between particles is of the order of the 
transition wavelength. Their effect is to disturb the equi- 
librium conditions of the cooling scheme, making it less 
effective than in the one-particle case and acting as an 
independent heat source with an intrinsic heating rate. 
Thus, trials with sideband cooling in optical lattices [5[ 
never reached sufficiently low temperatures to guaran- 
tee condensation. Our proposal is predicted, neverthe- 
less, to reach null temperature up to second order in the 
Lamb-Dicke parameter with low effect of the parametric 
fluctuation. In addition, this scheme has a faster cool- 
ing rate. This would facilitate counteracting the effect of 
dipole-dipole interaction by cooling faster and to a lower 
temperature, and make the goal of Bose Einstein conden- 
sation by laser cooling a step closer. 

Adapting the implementation proposal of FigE] to a 
trapping and cooling set up following [5] is relatively 
straightforward. Since the beams reproducing coupling B 
are highly detuned, they can also be used as the trapping 
light. Instead of a pair of counterpropagating beams, a 
set of three coplanar beams forming angles of 120^ with 
each other will generate an hexagonal dipole lattice. Ge- 
ometric considerations demonstrate that the projection 
of the three wavevectors on any axis of the lattice plane 
yields two wavevectors of opposite sign and same mod- 
ule, so that any axis on the trap plane will see the same 
pair of two counterpropagating beams of FigEl 

Although the harmonic wells forming the lattice are 
not strictly isotropic, they can be considered so in the 
neighborhood of the minimum. This is a correct assump- 
tion bearing in mind that the system resides in the Lamb- 
Dicke regime. This simplifies the implementation of the 
A coupling to the extent that it will have the same form 
of coupling B but for the fact that its plane needs to 
be tilted with respect to the trap in an angle such that 
satisfies condition in Eql6l Starting from a configuration 
parallel to that of coupling B, the tilting needs to be done 
with respect to the perpendicular of one of the trapping 
beams. 

This setup has the advantage that the axis perpendic- 
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FIG. 12. Numerical simulation of the steady state population 
for a group of 2 and 3 atoms. (Qs = z^, = 8.5z^, F = lOz^, 
A 2.5i/). 



ular to the trap will only be mechanically affected by the 
beams of the A coupling, so that sideband cooling can 
be implemented on that axis. The possibility of imple- 
menting the present cooling proposal also on this vertical 
axis would only require two counter-propagating beams 
of type B parallel to the vertical axis. 

This combination of cooling and trapping light con- 
verts it into a fairly simple experimental setting with a 
promising performance. 

A numerical simulation following the approach of [26| 
has been performed where the scheme was applied to a 
group of 2 and 3 atoms, obtaining results that differ lit- 
tle from the case of a single particle. Both dipole-dipole 
and superradiance terms were taken into account. The 
main difference with the single atom model is the pres- 
ence of the so-called bosonic-enhancement factors. The 
presence of many atoms in the same state enhances the 
couplings between levels. This has the effect of increas- 
ing the rate of interaction, but also of broadening ab- 
sorption and emission lines. For schemes relying on res- 



onance effects, like sideband cooling, this broadening is 
very detrimental. In our case, both carrier transition and 
blue sideband are under the influence of dark states and 
hence are not affected by this broadening. This is shown 
in Figdll where no difference in the steady state occu- 
pation values shows for increasing number of atoms. For 
this, the experimental values appearing in [5] have been 
taken. Since the trapping light is the same as that of 
coupling B, Qb — ^- This requires tjb = 4?]/^, which 
corresponds to a tilting angle of 60^. Cesium atoms can 
be used for which the linewidth of the excited state is 
F = lOv. We believe this scheme could hence reduce the 
reachable temperature by laser means in atom clouds, 
bringing BEC by all-optical means a step closer. 
XI. CONCLUSIONS 



A detailed study of a robust and fast laser cooling 
scheme for trapped three level systems has been pre- 
sented. It has been shown that a particular setting of the 
Lamb-Dicke parameters ratio combines two underlying 
cooling schemes so that the eigenspace dimension of final 
steady states is reduced to one. Thus, the ground state is 
reached up to second order in the Lamb-Dicke expansion. 
Being both a fast and a robust cooling scheme, two dif- 
ferent experimental implementations have been proposed 
that equivalent ly perform the scheme. An equally simple 
implementation would be able to also cool the two resting 
axis, thus converting it into a 3D cooling scheme without 
further changes in the experimental setup. It has been 
demonstrated that the robustness and fast cooling rate 
would make it suitable for cooling of several ions in a 
trap or of a cloud of atoms in an optical lattice, bringing 
BEC condensation with laser cooling closer to a factible 
goal. 
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Appendix A: Diagonalization of two Discrete States 
and a Continuum 



Let us first consider the manifold 
{|e) , {|+, 1)}, {| — , 1)}}. This can be diagonalized 
to generate the continuum Through |e), the 

bright state |+) also couples to the continuum 

with the term ^ (e|/c') a^^'^ ^ Let us now consider 
how the excited state overlaps with the continuum 
{\k)} originated in the manifold The 
Hamiltonian Hk that appears diagonal to \k) can be 
split into Hk'^ which is diagonal for \k') and the coupling 
to the state |+). It satisfies: 



Hk\k) 



k\k). 



\k) 



(Al) 



Prom ()A2p the following expression is obtained: 

1 



{k'\k) 



^/2k-k' 

which can be substituted in (|A3p to obtain: 



|(e|fc')l 



^-dk' = (fc-l^s). 



k-k' 

Now the overlap {e\k) can be calculated: 

(e|fc) = j {e\k') {k'\k)dk' 

= :^{k-nB){+\k) 



(A4) 



(A5) 



(A6) 



Appendix B: Effective Coupling of the Ground 
States 



In this section the relationship between the physical 
Lamb-Dicke parameter rjp and the effective t/b is derived. 
The method in |27| is followed, where the highly de- 
tuned excited level is adiabatically eliminated. Our tar- 
get Hamiltonian describing the coupling of both ground 
levels is the following: 



Vssh = ^{a(JM+VBa^;'^^) (Bl) 

The experimental implementation involving a highly 
detuned Raman coupling is described by the following 
interaction picture Hamiltonian: 



Hi = ^(|e>(t |e'^*(l + iVpib^e-'-^* + be"-')) 

+ \e){i |e'^*(l - irjpibh-'''' + be'"*)) + h.c). 

(B2) 

This Hamiltonian consists of the following harmonic 
terms: 

/^i = ^(|e){t|-|e>U|)6t 



/^2 = |^(|e>(t| + |e)a|) 
/^3 = ^(|e){t|-|e>U|)6 



(B3) 



Projection onto \k') and |+) yields: 



k'{k' 



V2 



{k'\e){+\k) = k{k'\k) (A2) 



Qb {+\k) + ^ / {k'\k) dk' = k {- 



(A3) 



with frequency values uji = A — uj2 = A and ujs = 
A-\-u. The derivation of the effective Hamiltonian follows 
the formula: 

3 

Heff{t)= ^ ow"^ [/i];,,/in] exp[i(a;^-a;n)^], (B4) 

m,n=l 
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where Umn is the harmonic average ujnin~^ = 
0J~^). Of the nine terms originating from this expres- 
sion, four of them are of second order in the Lamb-Dicke 
parameter. The remaining six can be expressed as: 



1 



+ 



4 A(A2 - iy2) 



(B5) 



+ 



4 A2 



(t4) 



where p = i{h — b^). This Hamiltonian is equivalent to 
()B1|I but for the atomic Stark shifts and the last term, 
which is proportional to l/A^. 

As long as A 3> i^, iJe/ / contains the target interaction 



Vssh with the parametric conditions: 



2 

0.3118 



p 

4A 

^2 2A 

Ja"^" A2 - i.2 



(B6) 



Solving for the Lamb-Dicke parameter: 



Vb = lip- 



2A2 



A2 



(B7) 



which yields rjB = ^rjp for the assumed hmit /S. ^ v. 
If we take the fluctuations optimal point = the 
condition becomes rjp/rjA = 2. This can be achieved for 
a layout where beam p is colinear to the trap axis and 
beam A, with a very similar frequency, is 60^ away from 
the axis. 



